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. Abstract. This paper is the first in a series of articles devoted to the construction and analysis of three-component 

' vector velocity fields in the gaseous disks of spiral galaxies, and to the discovery of giant anticyclones near 

, corotation which were predicted earlier. We analyse the line-of-sight velocity field of ionized gas in the spiral 

■ galaxy NGC 157 which has been obtained in the Ha emission line at the 6m telescope of SAO RAS. The field 

' contains more than 11,000 velocity estimates. The existence of systematic deviations of the observed gas velocities 

■««..^ , from pure circular motion is shown. A detailed investigation of these deviations is undertaken by applying a 

• recently formulated method based on Fourier analysis of the azimuthal distributions of the line-of-sight velocities 
at different distances from the galactic center. To restore the three-component vector velocity field from the 

Q . observed line-of-sight velocity field, two assumptions were made: 1) the perturbed surface density and residual 

' velocity components can be approximated to the formulae d cos(299 -I- Fi), where d and Ft are respectively the 

I amplitude and the phase; 2) the perturbed surface density and residual velocity components satisfy the Euler 

. equations. The correctness of both assumptions is proven on the basis of the observational data. As a result of the 

' analysis, all the main parameters of the wave spiral pattern are determined: the corotation radius, the amplitudes 

and phases of the gas velocity perturbations at different radii, as well as the velocity of circular rotation of the 
disk corrected for the influence of the velocity perturbations connected with the spiral arms. Finally, a restoration 
' of the vector velocity field in the gaseous disk of NGC 157 is performed. At a high confidence level, the presence 

of the two giant anticyclones in the reference frame rotating with the spiral pattern is shown; their sizes and the 
localization of their centers are consistent with the results of the analytic theory and of numerical simulations. 
Besides the anticyclones, the existence of cyclones in any residual velocity field is predicted. In the reference frame 
of the spiral pattern, detection of such cyclones is possible in the galaxies with a radial gradient of azimuthal 
residual velocity steeper than that of the rotation velocity. 

Key words, methods: data analysis - galaxies: individual: NGC 157 - galaxies: ISM - galaxies: kinematics and 
dynamics - galaxies: spiral - galaxies: structure 



1. Introduction 

The present work is devoted to the analysis of line-of-sight 
velocities of the gas in the spiral galaxy NGC 157. We at- 
tempt to evaluate all three vector components of the per- 
turbed gas velocity in the spiral density wave to find the 
angular velocity of the spiral pattern and the resonance 
locations, to obtain the rotation curve corrected for non- 
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circular motions connected with the wave structures, and 
also to demonstrate the existence of a vortex structure. 

Evidently, any method of restoration of the three- 
component vector velocity field from the observed one- 
component line-of-sight velocity field needs some addi- 
tional assumptions, and the correctness of the restoration 
completely depends on the adequacy of those assumptions. 
Here two assumptions were made: 1) the perturbed surface 
density and residual velocity components in the NGC 157 
can be approximated to the form d cos(2iy9 + Fi), where 
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Ci and Fi are respectively correspondingly the amplitude 
and the phase; 2) the perturbed surface density and resid- 
ual velocity components satisfy the Euler equations. In 
the paper we show that both assumptions agree with the 
observational data. First, the predominance of the second 
Fourier harmonic of the brightness over others allows to 
adequately present the perturbed surface density in the 
form a — Ca{r, t) cos[2ip — Fa{r)], whereas the prevalence 
of the second, third and the sine component of the first 
Fourier harmonics in the residual velocity field means that 
the perturbed velocity components were chosen correctly, 
in a similar manner. Second, analysis of the radial de- 
pendencies of the harmonic phases gives evidence for the 
existence of a connection between the line-of-sight veloc- 
ity disturbances of the gas and the observed spiral arms, 
in agreement with what follows from the Euler equations. 

In Sect. 2, a description of the observations of the 
galaxy is given. In Sect. 3, the results of an analysis of 
the line-of-sight velocity field under the approximation of 
pure circular motion are presented; we compare these re- 
sults with the earlier published observational data and 
discuss residual line-of-sight velocities obtained by sub- 
tracting the model velocity field from the observed one. 
We find that the deviations from pure circular motion 
demonstrate a systematic behaviour. In Sect. 4, a sim- 
ple model of the line-of-sight velocities is described, which 
takes into account motions in the two-armed density wave. 
It is shown that in this case the first three Fourier harmon- 
ics dominate throughout the line-of-sight velocity field. In 
Sect. 5, we present the results of the harmonic analysis of 
the observed line-of-sight velocity field which prove that 
the perturbed surface density and residual velocity com- 
ponents satisfy the Euler equations (i.e. the proof of the 
wave nature of the spiral structure). After that, restora- 
tion and analysis of the vector velocity field of the gas is 
performed. In Sect. 6, we show that the Fourier analysis 
of the line-of-sight velocity field allows us to restore the 
velocity field in the plane of the galaxy with a high degree 
of confidence. In Sect. 7, the corotation position is deter- 
mined directly from the observations. It is shown that the 
velocity field in the reference frame rotating with the spiral 
pattern clearly demonstrates the presence of two banana- 
like anticyclones, the centers of which lie near the corota- 
tion. The last Section contains the main conclusions. 



2. Observations and primary data reduction 

NGC 157 (Fig. |l|) is a rather bright {Bt = 11.00) and 
close (D ^ 21 Mpc for Hq = 75 km s"^ Mpc"^) isolated 
galaxy of SAB(rs)bc type (de Vaucouleurs et al. 1991). Its 
well-developed Grand-Design spiral structure is classified 
as Arm class 12 by Elmegreen & Elmegreen (1984). This 
galaxy is characterised by mildly enhanced star forma- 
tion, which follows from its rather high brightness in the 
Ha line and from far infrared flux densities (Ryder et al. 
1998). The large angular size (the isophotal radius Do/2 
= 2. '13, or ~ 13 kpc) and the moderate disk inclination 




-60 -40 -aO 20 40 60 

X offset (arcsec) 



Fig. 1. A monochromatic image of NGC 157 in the Ha 
emission line. North is roughly to the left (PAj^ — —76°). 
The phase curve of the second Fourier harmonic of the Ha 
brightness map is superimposed. Asterisks show the az- 
imuth positions of the maxima of the second harmonic at 
each radius. Positions of the maxima of the second Fourier 
harmonic of the NGC 157 if-band brightness map (kindly 
provided by S.D. Ryder) are also shown by squares. 

(axis ratio b/a — 0.65) make this galaxy very convenient 
for kinematic and photometric measurements. 

The optical morphology of this galaxy is rather com- 
plex. A blue image (Lynds 1974) shows a dusty and more 
or less symmetrical disk with two principal spiral arms 
and without any hint of a bar. Nevertheless there is some 
evidence of a weak bar, lacking bright HII regions, in the 
inner part of the disk, noticeable in the red (i?-band) im- 
age of the galaxy (Sempere & Rozas 1997). Note, however, 
that this bar does not reveal itself in the kinematics of Hi 
(Ryder et al. 1998). The spiral structure of NGC 157 is 
rather regular up to about 1', or ~ 6.1 kpc, from the cen- 
ter, and, by sight, becomes flocculent or multiarmed at the 
periphery. As follows from the Ha image of the galaxy, the 
two principal arms are not symmetrical with respect to the 
center (Sempere & Rozas 1997). According to Elmegreen 
et al. (1992) such an asymmetry may be responsible for 
the additional m=l component of the spiral wave struc- 
ture, driven by the two-armed spiral in the inner disk. 
Note also that one of the two predominant spiral arms is 
brighter in Ha than the other. 

Radio observations in the Hi line and at the 1.4 GHz 
continuum (Ryder et al. 1998) show that there is a close, 
though not strict correspondence of radio brightness and 
Ha surface brightness. In addition, NGC 157 appears to 
possess an extended warped Hi disk, which stretches out 
well beyond the optical borders. The rotation velocity of 
the Hi disk amounts to about half the maximal velocity 
of the optical disk, which enables us to conclude that the 
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galaxy has (1) a falling rotation curve and (2) an extended 
low-density halo beyond r ~ 8 — 10 kpc . The latter con- 
tains a small fraction of the total mass (about 10%) within 
the isophotal radius, so its gravitational input in the re- 
gion of the regular spiral structure is quite negligible. 

Rather bright Ha emission is observed not only inside 
the spiral arms, but also between them. This circumstance 
is favourable for the kinematic study of the ionized gas. 
The first investigation of the disk kinematics was carried 
out by Burbidge et al. (1961) who derived a slowly rising 
rotation curve from two long-slit spectra. Later Zasov and 
Kyazumov (1981), using long-slit spectral observations at 
different position angles, found a significant drop of the 
rotational velocity at r ~ 50"-55" in the north-eastern 
part of the disk. Recent Hi observations (Ryder et al. 1998) 
have confirmed the reality of this peculiar motion and have 
given evidence that it is caused both by a real declinion of 
the rotation curve and by asymmetry in the line-of-sight 
velocity distribution along the given positional angle. The 
optical rotation curve of the inner part of this galaxy was 
also obtained by Afanasiev et al. (1988). It was found that 
the rotation curve slowly rises, with Kot ~ 100 km s~^ at 
r = 1 — 2 kpc and reaching a fiat maximum of ~ 200 km 
s~^ at r ^ 5 — 6 kpc. 

The well-defined spiral structure in the inner angular 
radius of one arcminute makes NGC 157 suitable for the 
detailed investigation of density wave propagation in the 
disk. Analyzing the positions of different morphological 
features in this galaxy, Elmegreen et al. (1992) found that 
corotation takes place at about r ^ 56", so that the coro- 
tation circle limits the region of regular spiral structure. 
Note, however, that the determination of the location of 
resonances by morphological tracers involves a high de- 
gree of uncertainty (see the discussion in Sempere & Rozas 
1997). The alternative approach used by Sempere & Rozas 
(1997) is based on numerical simulations of the motions of 
interstellar molecular clouds in the stellar disk potential, 
as calculated from the i?-band image. The resulting distri- 
bution of molecular clouds has then been compared with 
the observed spiral structure of the galaxy. Although the 
results are model dependent, the obtained pattern speed 
in the best fit model is Vlp =40 km s^^ kpc^^, which cor- 
responds to a corotation radius of 50", close to the value 
claimed earlier by Elmegreen et al. (1992). 

Here we use a quite different method, the most direct 
one, to investigate the interconnection between the spiral 
wave and the kinematic behaviour of the interstellar gas. 
Our method is based on a Fourier analysis of the observa- 
tions of the line-of sight velocity field (Lyakhovich et al. 
1997, referred hereafter as L97). Some preliminary results 
of this work were given in Fridman et al. (1997, referred 
hereafter as F97). 

To obtain more detailed information about the kine- 
matics of gaseous clouds in NGC 157, we have undertaken 
observations in the Ha emission line with the scanning 
Fabry-Perot Interferometer (Dodonov et al. 1995) at the 
6-m Telescope of the Special Astrophysical Observatory 
of the Russian Academy of Sciences on October 24, 1995. 
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Fig. 2. The line-of-sight velocity field of the gas in 
NGC 157 with isovelocity contours overlaid. Systematic 
deviations from pure circular motion, connected with the 
spiral arms, are clearly seen. 

The Interferometer was mounted in a focal reducer (0"46 
per pixel scale) and installed at the prime focus of the 
telescope. To obtain 32 exposures 180 seconds each, we 
used the IPCS (512 x 512 pixels) as a detector (Afanasiev 
et al. 1986, 1995). Atmospheric conditions were pho- 
tometric with a seeing of 2.5 arcsec. The Fabry-Perot 
Interferometer was tuned at the 501st order and gave a 
velocity sampling of 19 km s~^. The order-separating fil- 
ter with bandpass FWHM of 10 A (A^ = 6613 A) was 
used to extract the chosen etalon order. Phase calibration 
was made using a neon line at A = 6598.95 A. 

The data were reduced with the ADHOC software 
package, developed in Marseille Observatory (Boulesteix 
1993). The resulting data cube of 256 x 256 x 32 pixels 
was created with a spatial sampling of 0.92 arcsec and a 
velocity sampling of 19 km s~^. The spatial resolution of 
2.5 arcsec corresponds to 250-300 pc, at the distance of 
NGC 157 with = 75 km s^^ Mpc"^ The line-of-sight 
velocity field as reduced from the data cube is shown in 
Fig. 1^. The isovelocity contours are superimposed on this 
image. Even on this figure one can clearly see the system- 
atic perturbation of the velocity field connected with the 
observed spiral arms. That becomes possible, in the first 
place, due to the continuous character of the velocity mea- 
surements, filling both the spiral arms and the interarm 
regions, which is what makes this particular galaxy very 
convenient for a detailed analysis of the velocity field. 

3. Pure circular motion model and residual 
velocities 

Systematic distortions of the line-of-sight velocity fields 
of a gas, induced by galactic spiral structure, are well 
known from observations (Lin et al. 1969; Yuan 1969; 
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Fig. 3. The rotation curve of NGC 157 obtained in the 
frame of the pure circular motion model for a set of fixed 
parameters of the galactic disk. Bars show errors at the 
level of 3 CT. 

Visser 1980). The amplitude of such distortions depends 
on a lot of physical and geometrical parameters, partic- 
ularly, on the wave amplitude, the form of the rotation 
curve, and the corotation position. The parameters listed 
above are to be determined. 

The usual approach to the problem was to deal sep- 
arately with the "unperturbed" rotation curve and the 
residual velocity field (obtained supposedly upon subtract- 
ing the said curve from the whole picture). But the prob- 
lem thus reformulated resisted all attempts to solve it in 
a direct way - in fact, it appears to be ill-conditioned. 
Indeed, to solve the problem as stated, one would have 
to know with high accuracy the unperturbed velocity at 
a given radius, which, in turn, is impossible without the 
knowledge of the perturbation velocities. The truth is, the 
problems of determining the rotation velocity and the per- 
turbations cannot be separated, because the perturbations 
caused by a density wave, being non-random, affect the 
shape of the restored rotation curve (see below Eqsjl], ^ 
and discussion). The strong interdependency of both parts 
of the problem is where the standard methods lose their 
ground (L97; F97). 

As the first step, we shall find the rotation velocity 
curve of the galaxy assuming the perturbations to be neg- 
ligible. It helps to determine the deviations from pure ro- 
tation, to see whether their behaviour is systematic or 
random, and in the former case to analyse the connection 
of the deviations with the spiral pattern. 

In the frame of this model the observed line-of-sight 
velocity can be written as (L97): 
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r) cos (p sini , 



(1) 



where v]/"iodi ^^le deviation of the observed velocity from 
the one predicted in the model of pure circular motion; Vs 
is a systemic line-of-sight velocity of the galaxy; 1^™°'^^ 
the model circular velocity of the gas in the plane of the 



galaxy which depends on the galactocentric radius r; (p is 
the galactocentric azimuthal angle measured with respect 
to the major axis; i is the inclination (the angle between 
the galactic plane and the sky plane) . The usage of Eqj^ 
gives the real rotation velocity of the galaxy if deviations 
yjymodi g^j.g random. Otherwise the value V^.^t"'^^ '^i^ differ 
systematically from the true equilibrium rotation velocity 
V,ll (L97, F97, see also EqsJ^, and Fig. |l| below). 

For the sake of definition, hereafter we identify the po- 
sition if — with the major semiaxis in the receding half 
of the galaxy and suppose the inclination to be negative 
if the angular rotation momentum is directed toward the 
observer. In this work we have assumed that the center 
coordinates, the inclination, and the line-of-nodes posi- 
tion angle are constant parameters which do not depend 
on the radius. Note that variations of the disk orientation 
parameters along the radius (often obtained as a result of 
the "splitting into rings" method of calculation of the said 
parameters, when each ring is considered independently), 
may be an artefact arising from neglecting the system- 
atic deviations of the gas velocities from the circular ones 
(L97). 

To determine the orientation parameters of the rota- 
tion plane of NGC 157 and its rotation curve, we have used 
more than 11 000 individual velocity measurement^ The 
model parameters for the galaxy were calculated by the 
method of least squares, minimizing the r.m.s. deviation of 
the observed velocities from the model ones derived from 
equation (|l|) . To estimate the rotation velocity at a radius 
r, we used the line-of-sight velocity values in the range 
r — d/2 to r -I- d/2. The value of the ring width d was cho- 
sen on the basis of the following arguments. Increasing d 
decreases the random errors, but at the same time it leads 
to the growth of the systematic error caused by neglecting 
gradients in the rotation velocity. The latter error could 
be estimated as Ssys ~ Kot d'^/^- Thus demanding the sys- 
tematic error to be less than 2 km s~^ and evaluating 1/"^ 
« (200 km s~^)/(20")^ we come to the restriction on d 
from above: d < 6" . On the other hand, it is meaningless 
to choose d smaller than the seeing value (2.5"). For these 
reasons, in the analysis presented below we use d = 5 pix- 
els or 4.6 arcsec. The rotation velocity was calculated at 
the step of one pixel (0.92 arcsec). After running through 
the whole range of the radial distances, we fit the rotation 
curve by a cubic spline. 

As a result, we have obtained the following best fit pa- 
rameters for NGC 157: the inclination i = —51.5° ± 5°, 
the position angle of the major axis (line of nodes) PAq — 
223.5° ± 2° (we identify PAq with the orientation (p = 
0), and the systemic heliocentric velocity of the galaxy 



^ It might be useful to emphasise here the difference between 
"spatially resolved" and "statistically independent" observa- 
tional points. In our case, the latter correspond to "individual 
velocity measurements", numbering up to 11 000. Evidently, 
the number of the "spatially resolved" observational points de- 
pends on seeing, while the number of the "statistically inde- 
pendent" observational points does not. 



A.M. Fridman et al.: Gas velocity field in NGC 157 



5 




-60 -40 -20 20 40 60 80 

X offset (arcsec) 



Fig. 4. The field of the residual line-of-sight velocities in 
NGC 157 with isovelocity contours overlaid. One can see 
that the deviations from the pure circular motion demon- 
strate a systematic character. 

Vsys = 1667 ± 5 km s^^. For comparison, in the RC3 (de 
Vaucouleurs et al. 1991) one can find the following esti- 
mates of the systemic velocity for NGC 157: 1730 ± 27 
km from optical data and 1668 ± 6 km s^^ from ra- 
dio (Hi line) data. Among the above listed numerical pa- 
rameters the value of i was estimated less reliably. From 
the isophote ellipticity the inclination is known to be 
i = 45° - 48° (Grosb0l 1985; BottineUi et al. 1984). In 
the paper of Ryder et al. (1998), where the same kinemat- 
ical data were used for the inner part of the galaxy and 
a more traditional method was applied to find the disk 
parameters, the inclination was found to be z = 45° ± 5°. 
Thus within the limits of observational errors our determi- 
nation of the inclination i agrees with the results of other 
authors. 

The obtained model rotation curve V^'^^'^^ is presented 
in Fig. ||. The errors of the rotation velocity (3 a) are 
shown by bars and are of order of 3 km s~^ (this value 
represents random errors only, systematic errors caused 
by the neglecting of motions in the density wave could be 
much higher). One can see that over a wide radial range 
(1 — 5 kpc) V™°'^^ rises quite linearly with the radius r; the 
fact was pointed out earlier by Barbridge et al. (1961) and 
Zasov & Kyazumov (1981). This peculiarity distinguishes 
NGC 157 from the majority of the other spiral galaxies, 
for which such regions of monotonous rise of the rotation 
curves rarely extend beyond 2 — 3 kpc. 

4. The model taking into account motions in the 
spiral arms 

In Fig. ^ we present a map of the gas residual veloci- 
ties (the difference between the observed velocities and 



the model ones obtained from (|l|)). As one can see, the 
observed velocities deviate systematically from the pure 
circular motion. The regions of maximum deviations have 
a spiral shape, which proves their connection with the spi- 
ral density wave in this galaxy. It means that the assump- 
tions underlying the model of pure circular motions are 
not valid and a more detailed model should be build. 

To construct a model describing the velocity field of 
NGC 157 more appropriately, it is necessary to take into 
account motions in the density wave. In the general case, 
the equation, connecting the components of the full veloc- 
ity vector in the galactocentric frame, Vr{r,ip), Vip{r,ip), 
and Vz (r, if) , with the observed line-of-sight velocity may 
be written as follows (L97; F97): 

V°^''{r, tp) = Vs + V^{r, ip) cos ip sini -|- 

-I- V;(r, ip) sin ip sini + Vz{r, ip) cost + (2) 

where ^1'°'^*' are errors in the observations. Values V°^^, Vg, 
Vr, and Vz describe the velocity field at the moment 
of the observations. In other words, they are functions of 
the galactocentric radius r and the azimuth ip but do not 
depend on the time. In the case of pure circular motion, 
V^r, = Kr'^'W, Vrir, ip) = 0, Vz{r, ip) = 0, so this 
formula evidently reduces to (|l|). 

The relation (||) is the only one between the line-of- 
sight velocity field and the components of the vector veloc- 
ity field of the gaseous disk, which can be derived without 
loss of generality. As a result, the direct restoration of the 
velocity field of the gaseous disk from the line-of-sight ve- 
locity field is impossible and some indirect methods should 
be apphed (F97). 

Particularly, one may pose the question in the follow- 
ing way. Let us take into account that we study a galaxy 
with clear two-armed spiral structure (for the general case 
of an m-armed spiral see L97; F97). As the arms are easy 
to trace, one may assume that the perturbed surface den- 
sity of the gaseous disk a (at the moment of observation) 
can be approximated by the expression 

a =. a(r) cos[2^-F,(r)], (3) 

where Co- and Fa- are the corresponding amplitude and 
phase. 

The validity of using (^ may be checked by Fourier 
analysis of the optical image of the galaxy: the approach 
is valid if the amplitude of the second Fourier harmonic 
predominates over the others and its phase curve closely 
correlates with the observed spiral. As one can see from 
Fig. H, the second harmonic dominates. The high level of 
the first harmonic in the Ha brightness map is caused by 
the nonsymmetrical distribution of star forming regions 
in the observed spiral arms. In our opinion, it does not 
reflect the real contribution of the first Fourier harmonic 
into the distribution of the perturbed surface density in 
the gaseous disk of the galaxy. The relatively low level 
of the first harmonics in the ii'-band image proves this 
assumption. Fig. |l| demonstrates that for NGC 157 the 
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Fig. 5. Averaged over the radial range 5" - 75" the contri- 
bution (dispersion given in arbitrary units) of the individ- 
ual Fourier harmonics into the deviation of the brightness 
distribution in different bands from the azimuthal symme- 
try, (a) Ha brightness map, (b) i^-band brightness map 
(kindly provided by S.D. Ryder). 



phase curve of the second Fourier harmonic of the Ha 
brightness map is in good accordance with the observed 
spiral pattern as well as with the second harmonic of K- 
band image of the galaxy. 

In a frame of the approach described, the perturbed 
velocity components may be written in a similar manner: 



Vrir, <p) = a(r-) cos[2ip - Fr{r)] , 
T4(r, <f)^C,{r) cos[2^-F,(r)]. 



(4) 
(5) 
(6) 



The validity of this approximation also can be checked 
directly from the observations. Indeed, taking into account 
that K-(r, ifl^ Vr{r, ^), F^(r, ^) = V.^^^^r) + V^(r, (^), 
Vzir, f) = Vz{r, (p), and substituting (|) - (|) in (||), we 
obtain the model representation of the line-of-sight veloc- 
ity (L97, F97): 



□ d2 



(r) cos3<f + b'^°'^^{r) sm3ip\ , 



(7) 



where Fourier coefficients related to the phases and am- 
plitudes of the velocity components are: 



mod2 



^mod2 



Cr sin Fr + Ctp cos F^ 



Cr COS Fr — Cip sin F„ 



mod2 

J mod2 
^2 

„niod2 



Cz COS Fj cot i , 
Cz sin Fz cot i , 
Cr sin Fr — 



Ccp cos Fcp 



"3 



d2 



Cr cos Fr 



Cp sin f y 



(8) 

(9) 

(10) 
(11) 
(12) 

(13) 



smi [af°'^'^{r) cos ip 



From the above relationships (g) - (jl^) we can see that 
the contributions of the different velocity components in 
the azimuthal Fourier harmonics of the observed line-of- 
sight velocity are distributed as follows. 

1. The systemic velocity of the galaxy contributes in the 
zeroth harmonic of the observed line-of-sight velocity. 

2. The circular rotation contributes in the coefficient of 
cosine of the first component of the observed line-of- 
sight velocity. 

3. The perturbed motion in the galactic plane contributes 
in the first and third harmonics of the observed line- 
of-sight velocity. 

4. The vertical motion (along the rotation axis) con- 
tributes in the second harmonic of the observed line- 
of-sight velocity^. 

Thus, if a galaxy has a pure two-armed structure, the 
observed line-of-sight velocity field should contain the ze- 
roth, the first, the second and the third harmonics only. 
Actually, if the structure is not pure two-armed as we 
have in the case of NGC 157 (see Sandage & Bedke 1994), 
and yet the two-armed mode dominates, one can expect a 
dominance of these harmonics in the observed line-of-sight 
velocity field (F97). 

^ Hereafter, by vertical motions we mean, first of all, the ver- 
tical motions in the density wave. The vertical velocity compo- 
nent in the density wave is antisymmetrical with respect to the 
central plane of the disk and thus leaves the disk flat on aver- 
age. These motions are observable due to the non-transparency 
of the disk in the Hq line. As a result, the contribution of the 
nearest part of the disk prevails. The connection of the ob- 
served vertical motions with the density wave can be checked 
analysing the correlation between phases of the second har- 
monic of the brightness map with the second harmonic of the 
line-of-sight velocity field (Sec. 5.3) or comparing amplitudes 
of the second harmonics of the line-of-sight velocity fields ob- 
tained in lines with different optical width (the method pro- 
posed in Fridman et al. 1998). 
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Fig. 6. Averaged over the radial range 5" - 75", the contri- 
bution of the individual Fourier harmonics into dispersion 
in the model of pure rotation. The prevalence of the first 
sine, second and third harmonics is clearly seen. 

Fig. |6| shows, for NGC 157, the contributions of the 
different Fourier harmonics to dispersion in the model of 
pure rotation averaged over the radius range 5" - 75". It 
is clearly seen that the deviations are mainly determined 
by the first (sine), second, and third harmonics. The indi- 
vidual contributions of any higher harmonic are relatively 
small. 

Hence, the observations of NGC 157 within the ra- 
dial range mentioned above justify the possibility to re- 
strict the input of motions induced by the density wave to 
the main {2ip) component only. It allows the use of rela- 
tions ^ - ([l3| ) connecting characteristics of the perturbed 
velocities and Fourier coefficients of line-of-sight velocity 
field to restore the three-dimensional velocity field of the 
galaxy. 

The results obtained in this way will be presented in 
Sec. ^. But before proceeding, we should check whether, 
indeed, there exists a close connection between the resid- 
ual velocities and the spiral arms observed in NGC 157. So, 
let us turn to the Fourier analysis of the line-of-sight ve- 
locity field to prove the wave nature of the spiral structure 
and the observed velocity deviations from pure rotation. 

5. Fourier analysis of the velocity field 

5.1. Fourier analysis technique 

Suppose that the inclination i, the position angle of the 
kinematical major axis (the line of nodes) PA, and the 
coordinates of the galactic center in the sky plane xq and 
j/o are known (say, having been determined in some inde- 
pendent way). 

If the observed line-of-sight velocity can be approxi- 
mated by a two-dimensional analytical function of galac- 
tocentric coordinates r and ip, then its expansion into a 
harmonic series may be presented as: 

y°bs ^Vs+smiJ2 {a°r!"{r) cos nip + bf^ir) sin nip){U) 

n=l 



where Umax is the number of the highest harmonic, to be 
chosen for each galaxy individually. 

A well-known technique has been applied to find 
Fourier coefficients at discrete points r = i?/ by minimis- 
ing the difference 

- J2 ian'"iRi)cosnipj+b'^'%Ri)sinnipj)smi\ .(15) 

n=l / 

To construct r. m. s. deviations x^(i?/) we use the ob- 
served values of Vj with the coordinates rj falling in the 
galactocentric radius range from Ri —d/2 to Ri+d/2 with 
d = 4.6 arcsec (see discussion on the choice of d in Sec. 3). 
Finally, to find a'^^lr) and b'^^{r) in intermediate points 
r we use the cubic spline drawn through the calculated 
points a°^''{Ri) and 6°''"(i?/). 

The recent investigation by Burlak et al. (2001) has 
shown that the method described here is very stable, and 
that its results are practically insensitive to the presence 
of 'holes' in the velocity fields if the filling factor of points 
with velocity measurements for a given galactocentric ring 
is greater than 30%. 

5.2. The large-scale line-of-sight velocity field 

The square root of the average dispersion in the model 
of pure rotation is about 19 km s~^, whereas the random 
error of a velocity measurement (in one pixel) is about 14 
km s~^. As the error (a) of the harmonic amplitude de- 
termination is about 1.5 km s~^, the above difference is 
statistically significant on the level of 99%. As mentioned 
above, deviation from the pure rotation is mainly deter- 
mined by the first (sine), second and third harmonics. 

The behaviour of the main harmonics with galactocen- 
tric radius is shown in Fig. Bars correspond to 3 cr level. 
This figure demonstrates that the amplitudes of the main 
harmonics exceed the errors substantially, and hence, their 
values can be determined from the observations reliably 
over most of the galactic disk. The first plot shows the ab- 
solute value of the first harmonic cosine coefficient^. Due 
to the good filling of the galactic disk by the velocity mea- 
surements, the value of a°^^ is close to the value calculated 
with no account for the contribution of the other Fourier 
harmonics, that is by the reverse Fourier transformation^. 
A typical discrepancy does not exceed 1 km s~^, which is 
less than the errors of the af'^ determination. The maxi- 
mum difference (of the order of 5 km s~^) takes place near 
the galactic center where there are only a few measure- 
ments suitable for the Fourier analysis, and at the edge 

The proper value of the coefficient a^^^ is negative, which 
corresponds to counterclockwise galaxy rotation from the ob- 
server's point of view. 

* The latter value is equal to rotation velocity Kot°'^^ calcu- 
lated in the model of pure circular motions. 
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Fig. 7. Radial behaviour of the main harmonics of the 
line-of-sight velocity field for NGC 157. Bars show the 
observational errors (3cr). 

of the galaxy where "holes" exist in the observed veloc- 
ity field. It should be noted that the mentioned closeness 
is a direct consequence of the independence of different 
Fourier harmonics and bears no relation to the validity of 
the model of pure circular motions. According to Eq. 
the difference between ai and the equlibrium rotation ve- 
locity is systematicaly nonzero and about the amplitude 
of motions in the density wave (i.e about 10-20 km s^^, 
see Fig. ^ below) . 

Now we restrict the analysed field of the line-of-sight 
velocity of the galaxy only to the first, second, and third 
harmonics: 

V'^'^/smi = Vsys/sini + 

+ af'^r) cos ^ + bl'"%R)smip + 
af%r) cos2ip + 6f"(r)sin2(p + 

+ af'^r) cosSif + bf^r) sinSip. (16) 

The averaged velocity dispersion in the model is about 
15 km s~^. The restriction to the first three harmonics 
is reasonable from the goodness-of-fit point of view (see 
Bevington 1975) because it enables the elimination of sys- 
tematic deviations of the observed velocity field from those 
expected in the case of pure rotation of the gas in the plane 
of the disk. 
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Fig. 8. Deviations of the observed line-of-sight velocity 
field from that restricted to the first three Fourier har- 
monics. The chaotic small-scale structure of the field is 
clearly seen. It is caused by random errors and by small- 
scale motions in the star-formation regions. 



Fig. 1^ presents the residual velocity field for this model 
- the observed velocity field minus the model one, re- 
stricted to the first three Fourier harmonics according to 
(|l6|). As one can see, the residual velocities demonstrate a 
small-scale chaotic structure. The areas of maximum de- 
viations correlate with the intense star-formation regions 
seen on the Ha image of the galaxy (Fig. |l]). 

Therefore we can conclude that the first three harmon- 
ics describe the large-scale structure of the velocity field 
rather well, while the higher-order harmonics over the ma- 
jor part of the galaxy have smaller amplitudes and relate 
mostly to small-scale distortions in the brightest star for- 
mation regions and to random errors. Hence for further 
analysis of the large-scale structure of the velocity field 
we can restrict ourselves to the first three harmonics. 

In the frame of the three-harmonic model we have 
reevaluated the parameters of the orientation of the galac- 
tic rotation plane {i and PAq). Our calculations show that 
the parameters which minimize the r. m. s. deviations of 
the line-of-sight velocities from (^6|) coincide with those 
found in the frame of the pure circular model within the 
measurement errors. This is easy to explain, remembering 
that the harmonic coefficients oscillate along the radius 
(see Fig. |^), so, being averaged, their input is close to 
zero. 

Therefore, for a given galaxy and for a given rich statis- 
tical sample of measurements of the line-of-sight velocity, 
the model of pure circular motion allows to find parame- 
ters of the orientation of the disk with a good accuracy, 
even without taking into account the motions in the spiral 
arms. 
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Residuals of the observed velocity field 
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Fig. 9. Large-scale component of the full (top) and resid- Fig. 10. The symmetric (top) and antisymmetric (bot- 

ual (bottom) line-of-sight velocity fields. Only contribu- tom) parts of the large-scale residual line-of-sight velocity 

tions of the first, second and third Fourier harmonics are field, 
included. 



Fig. H (top) shows a model velocity field of NGC 157 in 
the sky plane, calculated by using the formula (|l|). Fig. | 
(bottom) shows a large-scale residual velocity field which 
represents the model velocity field upon subtraction of the 
first cosine harmonic. The picture of this residual velocity 
field is rather complex. 

To clarify the structure of the velocity field, in Fig. ^ 
we present a symmetric (top) and antisymmetric (bottom) 
parts of the residual velocity field separately. The symmet- 
ric part contains the contribution of the second harmonic, 
and the antisymmetric part - contributions of the third 
harmonic and of the first sine harmonic. 

Note that although Fig. |o| (bottom) resembles the 
residual velocity maps of Canzian (1993), use of Canzian's 
method for the determination of the corotation radius be- 
came possible only after separating the residual velocity 



field into symmetric and antisymmetric components which 
has not been done by Canzian. As for the original field of 
residual velocities (Fig. ^, bottom), its shape is infiuenced 
by the second Fourier harmonic and does not contain any 
switching from one-armed to three-armed spiral, that is, 
does not allow use of Canzian's method in principle. It is 
also worth mentioning that Canzian (1993) did not take 
the second harmonic into account at all. 

5.3. Proof of the wave nature of the grand design arms 

As mentioned above, the dominance of the first, second 
and third harmonics in the expansion of the line-of-sight 
velocity field of NGC 157 into Fourier series finds a natural 
explanation under the assumption of the wave nature of 
the grand design structure (see also F97). 
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Fig. 11. Radial behaviour of the amplitude (top) and the 
phase (bottom) of the second Fourier harmonic of the line- 
of-sight velocity. Bars show the observational errors (Scr). 

Let us obtain direct evidence for the latter on the base 
of the connection between the main Fourier harmonics of 
the line-of-sight velocity field and the spiral structure ob- 
served in NGC 157. To do this we consider in detail the 
behaviour of the second and the third Fourier harmon- 
ics. For further consideration it is convenient to define the 
amplitude C„ and the phase Fn for any harmonic by the 
following equations: 



C„ cos(F„) 
C„ sin(F„) = 6f 



obs 



(17) 
(18) 



Fig. |Tl| presents radial variations of the amplitude and 
phase for the second harmonic of the line-of-sight velocity, 
following (|r^) and (|l8|) for n = 2. The bottom plot shows 
the locus of the second harmonic maximum (i.e. where 
cos[2(/3 — F2{R)] = 1) on the galactic plane, which demon- 
strates the behaviour of the second harmonic phase. Arcs 
on the plot show errors of the phase determination (3 cr). 
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Fig. 12. Superposition of the second harmonic of the line- 
of-sight velocity field on the Ha image of NGC 157 (gray 
scale). Triangles show the azimuth positions of the max- 
ima of the harmonic at each radius. The good correspon- 
dence of the phase curve with the observed position of 
the spiral arms gives strong evidence for the non-circular 
velocities, associated with the spiral structure. 

As one can see, there are only a few radial positions (near 
r = 20" and 50") where the amplitude of the second har- 
monic is comparable with the error of its measurement 
(C2 < 3 • (Tea), while over most of the galaxy it is deter- 
mined quite safely (C2 > 3 • (Tc2 ) • The lines of the second 
harmonic maxima look like a two-armed trailing spiral, 
and they are well correlated with the galactic spiral arms 
observed in the emission line of Ha inside 50" (Fig. |l2|) . 
This gives strong evidence that the second harmonic of 
the line-of-sight velocity is related to the vertical motions 
in the observed spiral density wave in the region r < 50". 
In the outer region, the second harmonic of the line-of- 
sight velocity can be caused by the plane motions in m = 
3 mode, which becomes relatively stronger in this region. 
A discussion of this subject is beyond the scope of the 
current paper. 

Fig. [ij presents the behaviour of the amplitude and 
phase for the third Fourier harmonic of the line-of-sight 



velocity field, which are defined by equations (17) and 
(|l8|) for n = 3, as well as their errors. The bottom plot 
shows the locus which is determined by the equation 
cos[3((5 — F3{R)] = 1 on the galactic plane. The amplitude 
of the third harmonic over most of the galaxy is deter- 
mined rather reliably (C3 > 3 • ctcs)- The lines of maxima 
of the third harmonic also look like trailing spirals over 
most of the galaxy, excluding a small area at the periph- 
ery (r > 60"), where they have the form of the leading 
spirals. 

Though the spiral-like behaviour of the third harmonic 
is evident, a direct comparison of its pattern with the ob- 
served spiral arms of the galaxy is impossible because of 
the unequal number of arms. Nevertheless, it can be shown 
that if the observed spiral structure is a tightly wound den- 
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Fig. 13. Radial behaviour of the amphtude (top) and the 
phase (bottom) of the third Fourier harmonic of the hne- 
of-sight velocity. Bars show the observational errors (Scr). 

sity wave, there must exist a connection between the phase 
of the third Fourier harmonic of the linc-of-sight velocity 
field and the phase of the perturbed surface density (L97; 
see also the Appendix below): 

^ + 7r/2. (19) 

Within the corotation radius, this relation will be valid if 
C,p > Cr- Under the opposite condition C^p < Cr, F^ « 
Fa- — Tr/2 within the corotation radius. 

To check the validity of relation (p^, it is convenient 
to compare the locus of maxima of the observed surface 
density of gaseous galactic disk and two-armed spiral with 
the phase value equal to the phase of the third harmonic 
shifted by 7r/2. In Fig. ^ one can see the superposi- 
tion of the maxima of the modified third harmonic of the 
line-of-sight velocity field on the maxima of the second 
Fourier harmonic of the K-hand surface brightness map. 
The modified third harmonic varies as cos(2(y9-|-7r/2 — F3), 
where i^3(i?) is the phase of the original third harmonic. 
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Fig. 14. Superposition of the modified third harmonic 
of the line-of-sight velocity field (triangles) on the sec- 
ond Fourier harmonic of the K-hand brightness map of 
NGC 157 (squares). The modified third harmonic has a 
form of cos{2ip + 7r/2 — F^), where F3 is the phase of the 
original third harmonic. Points show the azimuth positions 
of the maxima of the relevant harmonic at each radius. A 
good agreement of the two phase curves gives strong ev- 
idence for the non-circular velocities, associated with the 
spiral structure. 

The good correspondence of the phase curve with the ob- 
served position of the spiral arms (see also Fig. ^ proves 
that the observed deviations of the line-of-sight velocity 
field from the pure circular motion on one hand and the 
spiral arms seen in the brightness map of the galaxy on 
the other hand are two different manifestations of the same 
phenomenon - the spiral density wave. 

It is worth noting that the transformation of the trail- 
ing spiral into the leading one at the periphery of the disk 
which is seen in Fig. |lj may be a real feature of the spiral 
pattern geometry, since this peculiarity appears both in 
the Ha and K images (see Fig. |]), as well as in the be- 
haviour of the third Fourier harmonic of the line-of-sight 
velocity field (Fig. 

6. Determination of the vector velocity field in 
the plane of the galaxy 

It was shown above that the large-scale structure of the 
velocity field may be well described in the frame of a 
model which includes the first, the second and the third 
azimuthal Fourier harmonics. It was also demonstrated 
that this result, as well as the structure of the line-of-sight 
velocity field, may be naturally explained if one suggests 
the residual velocities to be caused by the density wave 
dominated by a two-armed mode. These results allow the 
restoration of both the equilibrium rotation curve of the 
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galaxy and the velocity components in the density wave, 
by equating correspondent Fourier coefficients in Eqs.(^ 
and (|l|): 



Amplitude A^. 



^mod2 _L sin + COS 



+ 



^obs 



- acos^V + C^sini^^ = 2&f% 
C^cosi^^coti = af% 
C;,sinF^coti = 6f% 

- asinF^ + C^cosF^ = ^af"", 



Cr COS Fr 



C^sinF^ = 2 6f% 



(20) 

(21) 
(22) 
(23) 
(24) 
(25) 




R (arcsec) 

Amplitude B 



An additional difficulty still exists, since both rotation 
and motions in spiral arms contribute to the first cosine 
harmonic of the line-of-sight velocity. As a result of this in- 
terference, the system of equations is incomplete. Indeed, 
to determine seven unknown functions on the left-hand 
side of Eqs. (^ - (^ we have only six measured Fourier 
coefficients. 

From Eqs. (^) - ( ^5| ) it follows that the determina- 
tion of the coefficients of the Fourier harmonics of the 
observed line-of-sight velocity field gives a possibility to 
determine the parameters = Cz cos Fz, = Cz sin Fz , 
Ar = Cr COS Fr, and = C^sinF^, without any addi- 
tional assumptions. The first pair of these coefficients al- 
lows to restore the vertical velocity. Ar characterises the 
behaviour of the perturbed radial velocity along the dy- 
namical major axis of the galaxy, and B^ characterises 
the amplitude of the perturbed azimuthal velocity for the 
points at the minor axis. The radial dependencies of the 
coefficients Ar(r) and B^p{r) are shown in Fig. 

However the full set of unknown parameters: V^°'^^ 
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Fig. 15. Radial behaviour of the amplitudes of the co- 
sine component for the radial residual velocity (top) and 
of the sine component for the azimuthal residual velocity 
(bottom). 

The errors of determination of are the most cru- 
cial for application of the above relations. To avoid non- 
reliable estimates we assume the following condition for 

> 3 • fj^oba , where 



63 to be suitable for calculations: \h\ 



pbs 
3 



C^, F^, Cr, and F,. cannot be derived unambiguously ^bf- is an error of the measurements, 
without additional conditions. To close the set of equa- 
tions (|^) - (|2|), we use two independent approaches pro- 
posed in L97 and F97. 



Fig. 16 shows the rotation velocity (V^ 



mod2\ 



6.1. Method I, based on the relation between the 
phases of the radial and azimuthal residual 
velocities 

The first way to solve Eqs.(^-^) is to use the theoretical 
relation between the phases of the radial and tangential 
residual velocities in a tightly wound density wave (L97; 
F97): 



the ampli- 
tudes and the phases of the radial and tangential residual 
velocities calculated from the above equations in the re- 
gions where \h'^^\ > 3 ■ (Jjobs . For comparison, the rotation 

curve 1^™°'^^ obtained in the model of pure circular mo- 
tion (it corresponds to the Fourier coefficient af^^) is also 
plotted, as a dashed Hue. One can see from this figure 
that the maxima of the radial and tangential residual ve- 
locities outline the shape of trailing spirals, whose pitch 
angle is close to that of the visible spiral arms. Within 
7" < r < 43", the phase difference between the radial 
and tangential velocities Fr — F^ is equal to 90° , and out- 



F — F — 

^ r 



7r/2 at r < Tc 
— 7r/2 at r > Vc 



side 



43" Fr - F„ 



-90°. Therefore, this method 



sgn{Lj) Tr/2 . 



(26) 



By using this relation one can obtain from (^ - ^5|) 
(L97; F97): 



mod 2 



^obs 



Lobs „obs /Lobs 



A^ EE C^cosF^ = af"' {bf^b':^' 



Br 



Cr sin Fr 



-.obs fuohs /uobs 



1), 



(6f76§ 



1). 



(27) 
(28) 



enables us to localize the corotation radius in the region 
around r = 43" (see Eq. 26). The mean amplitudes of the 
perturbed radial and tangential velocities are 20 - 30 km 
s~^. The rotation velocity found here differs from the ro- 
tation velocity determined in the model of pure circular 
motion by amounts comparable to the velocity amplitude 
in the spiral arms. 
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Fig. 16. The rotation curve (dashed hne shows radial behaviour of af^^ which corresponds to the rotation 

velocity in the model of the pure circular motion, V^'^^'^^), the amplitudes of the radial and azimuthal residual velocities 
(Cr and Cip), the phase difference Fr — F^p, and the localisation of the maxima of the residual radial and azimuthal 
velocities at each radius, obtained by using the relation between the phases of radial and tangential residual velocities. 
The errors are shown at the level of 3a. 



6.2. Method II, based on the relation between the 
phases of the residual radial velocity and 
perturbed surface density 

Another way to close the system of Eqs.(^- ^5]) is to use 
the relation between the phases of the perturbed surface 
density and the radial residual velocity, in the approxima- 
tion of tightly wound spirals (L97; F97): 

^.-F. [i-sgnm./2. (29) 

The determination of the perturbed density phase us- 
ing the image of a galaxy is a separate problem. Here we 
assume that the locations of the maxima of the perturbed 
surface density at a given radius coincide with the loca- 
tions of the maxima of the surface brightness related to 
the spiral arms. This assumption allows to determine the 
behaviour of the phase of the second Fourier harmonic of 
the perturbed surface density Fa- by applying the Fourier 
analysis to the brightness distribution in a galaxy. The 
values of F^ and cos F^ obtained in this way from the K- 
band image of NGC 157 (see Fig. |) are shown in Fig. |l^. 

From (|^) - (HI) and (H) we have (L97; F97): 

y^^°d2 = flObs „ gObs _ (6°b- - 6f ^)tanF<, , (30) 



= 2af'^ + (6f ^ - 6f ^) tanF, , (31) 

Br = {bf' - 6f")tani^<,. (32) 

Fig. |l8| presents the radial dependence of the same pa- 
rameters as in Fig. |l^, but obtained by Method II for 
those regions where cos-Fo- can be determined most reli- 
ably, that is, where jcosF^I > 3 • (Tcosf„ (see Fig. |l^. As 
one can see from Fig. ^ the lines of the maxima of the 
radial and tangential velocities look like trailing spirals 
whose appearance agrees rather well with the phases ob- 
tained earlier by Method I (see Fig. |l^). The change of 
the phase difference between the radial and tangential ve- 
locities from Fr - F^ ~ 90° to Fr - F^ -90° occurs 
in the radial range 36" < r < 63". Therefore, the corota- 
tion radius obtained by this method is located within this 
range - in agreement with the estimate found above by 
Method I. The characteristic amplitudes of the radial and 
tangential velocities are of the order of 20 - 30 km s~^, 
again in good agreement with the previous estimates. 

6.3. Final model of the velocity field 

The most significant systematic error in the estimates ob- 
tained above arises from the assumption of a sharp switch 
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Fig. 18. AS Fig. 16, now obtained by using the relation between phases of the radial residual velocity and the perturbed 
surface density. The errors are shown at the level of 3a. 



of the phase relations near corotation. In reality, the area 
of switch cannot be too narrow. Our estimates show (L97) 
that the switch area size is of order of 1/kr, which for 
NGC 157 corresponds to 10"- 20". Therefore, the veloc- 
ity field parameters derived above have to be smoothed 
near corotation by using a filter with a width of about 

1/kr. 

Another argument in favour of the necessity of smooth- 
ing is the behaviour of the rotation curve near the coro- 
tation found by Methods I and II. Its non-monotonic fea- 
tures, seen in Figs. |l^ and 18, directly result from the 



assumed sharp switch of the phases, since the radial distri- 
bution of the azimuthally-averaged stellar disk brightness, 
which follows the surface density distribution, is rather 
smooth and hence gives no reason to expect any peculiar- 
ities of the rotation curve near r = 30" — 50". 

So we claim the necessity to smooth the calculated ro- 
tation curve near the corotation. That, in turn, leads to 
a broadening of the switch area where the phase differ- 
ences of the radial and tangential residual velocities and 
the perturbed surface density change. 

To obtain the best model of the velocity field for the 
galaxy, the radial dependence of the residual velocities was 
varied within the range of uncertainties given by the com- 
parison of Methods I and II. The "best fit" model was 
chosen as a model which gives a smooth shape of the ro- 



tation curve, which can be expected for a smooth mass 
distribution in the galaxy. The rotation curve for the best 
fit model is shown in Fig. |l^. For comparison, the rota- 
tion curve obtained in the traditional way, using a model 
of pure rotation, is also presented. In the latter case, the 
presence of wave-shaped features is seen. Similar wave-like 
details of rotation curves were found in many galaxies, and 
they can naturally be explained in the frame of the density 
wave theory (Lin et al. 1969; Yuan 1969). 

The estimation of the rotation velocity 1^.™°'^^ at dif- 
ferent radii enables us to determine unambiguously the 
other characteristics of the galactic velocity field from the 
set of equations (1 



20 shows the behaviour of 



_-|2^). Fig. 

the rotation curve, as well as that of the amplitudes and 
phases of radial and tangential velocities in the density 
wave, as obtained for the best-fit model. 



7. Corotation radius of the spiral pattern 

The principal aim of the present work is discovery of new 
galactic structures - giant anticyclones - near the corota- 
tion. To make these visible, knowledge of the corotation 
radius of the density wave is needed. Below we applied the 
methods proposed in L97 and F97 to determine the posi- 
tion of corotation from the observed line-of-sight velocity 
field of the ionized gas in the galaxy. 
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Fig. 20. The rotation curve V^i-ot"'^^ (dashed hne shows the radial behaviour of af'^, which corresponds to the rotation 
velocity in the model of pure circular motion, V^.™""^^), the amplitudes, the phase difference F,. — F^p, and the positions 
of the maxima of radial and tangential residual velocities at different radii obtained in the best-fit model (see the text). 



7.1. Method 1, using the relation between the phases 
of radial and tangential residual velocities 

From Eqs. (|^) - (|2|) and ^ 
be obtained (L97; F97): 



|26|) the following relation can 



Lobs 



uohs 
^3 



(6f ^ -r u. 

Hereafter, the upper sign in plus-minus and minus plus 
signs corresponds to the region inside the corotation, and 
the lower sign - outside it. Taking into account that the 
amplitudes Cr and have always positive values, it fol- 
lows from equation (|33|) that 



5obs)2 



(33) 



|6f"(r)| - |6f'^(r)| < 0, 
\bf^{r)\-\bf'^{r)\>0, 



for 
for 



r <rc, 

r > Tr. 



(34) 



These inequalities offer a possibility to determine 
the location of corotation from the observational data. 
According to Eq.(|3^), corotation is located in the region 
where the difference of the amplitudes of the third sine 
and the first sine harmonics (|63'^''(?')| — Wi^^ir)]) changes 
its sign from minus to plus. 

Fig. m shows the radial behaviour of 1 6!^'''' (r) | - 1 bf"" (r) | 
in NGC 157. As one can see from the figure, this function 
is negative within the errors in the inner part of the galaxy 
and positive in the outer region, in accordance with the 



expectations. From these data it follows that the corota- 
tion radius is about 42" ± 5". 



7.2. Method 2, using the relation between the phases 
of the radial residual velocity and the perturbed 
surface density. 



System (M) - (E5h closed by relation (gg) gives (L97; F97): 



obs 

'3 



LobsN 



cosF„ 



=F Cr cos F„ , 



from which follows bf that 
{bf'-ir) - 6f"(r)) cosF^(r) 



{bf%r) 



bf^{r)) cosF^(r) 



0, 
0, 



for r < Tc 
for r > Tr 



(35) 



(36) 



These two inequalities allow to find the location of corota- 
tion on the base of the brightness map and the line-of-sight 
velocity field. 

Fig. 22 shows the dependence of (63'^'' — bf^^) cos F„ 
on galactocentric radius r. The estimation of Tc by this 
method gives 42" ±6". 

It should be noted that the methods described above 
may possess not only random errors, but also some system- 
atic errors caused by the violation of the approximation 
of tightly wound spirals. However they use observational 
data related to different parameters of the spiral struc- 
ture, and their systematic errors must be different. So the 
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Fig. 17. Phase of the second harmonic of the perturbed 
surface density (assumed to be equal to the phase of the 
second harmonic of the K brightness map) and the value 
of cos Fa used to calculate the parameters of the velocity 
field according to Method II. 
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Fig. 19. Rotation curve of NGC 157 in the best fit model 
(solid line) and in the model of pure rotation (dotted line 
with error bars 3a). 
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Fig. 21. Trend of |6§''"(r)| 



3 y )\ ~ ^ith galactocen- 

tric radius r in NGC 157. Under the approximation of 
tightly wound spirals this difference should be negative 
inside corotation and positive outside it. According to the 
(lata presented, the corotation radius is about 42" ± 5". 
Error bars correspond to 3a level. The jumps of the ve- 
locity components in the regions r < 20" and r > 65" 
unveil some peculiarities of the galactic structure in these 
regions. In the central area it may be ascribed to the pres- 
ence of a bar, at the periphery these features are located 
where we observe a strong kink of spiral arms. 
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Fig. 22. Variation of {bf^^ — bf^'^) cosFa with galactocen- 
tric radius r in NGC 157. In the WKB-approximation the 
difference should be negative inside corotation and posi- 
tive outside it. Error bars correspond to 3a level. These 
data show the corotation radius to be about 42" ± 6". 



results of the application of these methods give evidence 
that the radius of corotation is within 36"-48"and that 
the characteristic value of the systematic error introduced 
by the tightly wound approximation hardly exceeds 6". 
For the second method described above, the error may 
be slightly higher because in this case an additional as- 
sumption about one-to-one correspondence between the 
azimuthal positions of the maxima of perturbed density 
and perturbed luminosity was made. Note that a slightly 
higher value for the corotation radius, r = 50", was re- 
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Fig. 23. Positions of the main resonances of the spiral 
structure in NGC 157 for the corotation position rc = 42". 



cently obtained from the numerical simulation of NGC 157 
by Sempere and Rozas (1997). 

Assuming the location of the corotation radius to be 
within the radius range of 36"-48", and taking the ro- 
tation curve as obtained in the best fit model, the inner 
Lindblad resonance (ILR) is either complitely absent in 
this galaxy, or located near the very galactic center, where 
the shape of the rotation curve is poorly determined. The 
outer Lindblad resonance (OLR) is located near the edge 
of the ionized-gas disk, in the region where the well-defined 
spiral arms end. Therefore, the bright spiral arms are con- 
strained to exist inside the OLR. This result agrees with 
the results of the numerical simulations of the spiral pat- 
tern in NGC 157 (Sempere & Rozas 1997). On the other 
hand, a strong kink of the spiral arms (change from trail- 
ing to leading spirals) is located near the outer 4:1 reso- 
nance (compare Fig. [l^ and Fig. 23). 

The estimate of the difference between the phases of 
the radial perturbed velocity and the azimuthal perturbed 
velocity in the best fit model allows to find the corotation 
radius more accurately by using the relation between the 
phases of radial and azimuthal velocity perturbations at 
the corotation. As was shown in L97, the phase difference 
between the radial and azimuthal perturbed velocities be- 
comes zero at corotation. In our case, this corresponds to 
Tc — 42". This value will be used in the next section to 
restore the velocity field of the gas in the reference frame 
rotating with the spiral pattern. 

It is worth noting that the corotation radius was found 
above by using some particular values of the disk inclina- 
tion angle i and the position angle of the line of nodes PAq, 
which give the minimum dispersion in the model taking 
into account the first, second, and third harmonics of the 
line-of-sight velocity field. Any changes of these parame- 
ters, even within the errors of their determinations, may 
affect the determination of the location of corotation. So 



we have carried out a special analysis of the influence of 
the variations of inclination and position angle on the es- 
timate of the corotation radius. We find that variations 
of the inclination by ±5° do not change (within the ac- 
curacy of 1"- 2") the position of the corotation radius. 
However the value of Vc strongly depends on the accepted 
PAo . Variations of PAq by ±2° result in corotation radius 
variations of about ^ 5". However, the qualitative picture 
of the restored velocity field remains the same even in this 
case. 



8. Vortex structures 

As a first step to visualising the restored velocity field of 
NGC 157, in Fig. ^ we show a vector field of the resid- 
ual velocities of the gas in the plane of the galactic disk, 
calculated for the best fit model. We can see four vortices 
near Vc - two cyclones and two anticyclones. This result is 
found for all values of parameters used for the restoration 
within the range of uncertainties. 

The restored velocity field in the reference frame rotat- 
ing with the pattern speed of the density wave, obtained 
in the best fit model, is shown in Fig. g5| The gas veloc- 
ities inside the spiral arms are directed preferably along 
the arms: toward the center inside the corotation circle 
and outward outside the corotation. It is clearly seen that 
the velocity field of the galaxy demonstrates two anticy- 
clones located in the corotation region between the spi- 
ral arms. Such anticyclonic vortices had been predicted 
on the basis of laboratory experiments on shallow water 
(Nezlin et al. 1986). The velocity amplitude of the vortices 
is about 30 km s~^. The two cyclones seen in the resid- 
ual velocity field are absent in the full velocity field in the 
reference frame corotating with the spiral arms. To make 
the cyclones visible, the radial gradient of the perturbed 
azimuthal velocity {dV^/dr) should be larger than the ra- 
dial gradient of the rotation velocity (d( Vrot — ^cf) / dr) in 
the same region of the disk. It is necessary to transform 
the anticyclonic shear, caused by the disk rotation, into a 
cyclonic shear. In the gaseous disk of NGC 157 this con- 
dition is not fulfilled. The cyclones, if they exist, must be 
confined to a small number of resolution elements in our 
data so their detection is very difficult. In galaxies with a 
stronger gradient of cyclones can be observed^. 

As shown by numerical simulations (Baev et al. 1987; 
Baev & Fridman 1988), spirals and vortices are gener- 
ated simultaneously by one and the same instability, al- 
ready at the linear stage. The levels of saturation of the 
instability are different for different perturbed functions 
- the densities of gaseous and stellar disks and their ve- 
locities. The amplitudes of the velocities and stellar disk 
density perturbations stopped growing at a linear stage, 
while the amplitude of the gaseous disk density perturba- 
tion becomes nonlinear. As has been shown many times 



Within the three years of the present paper publication the 
prediction was confirmed (Fridman et al. 1999; Fridman et al. 
2001a; Fridman et al. 2001b). 
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Fig. 24. Restored vector field of the residual velocity in gaseous disk of NGC 157. Overlaid squares for each galacto- 
centric radius show azimuthal positions of maximum value of the second Fourier harmonics of the Ha brightness map 
of the galaxy. 



during the past 15 years, the streamlines of liquid parti- 
cles in a self-consistent field of hydrodynamical and/or 
self-gravitational forces form cyclones and antycyclones 
(Nezlin et al. 1986; Baev & Fridman 1988; Afanasiev & 
Fridman 1993; Lyakhovich et al. 1996; F97). The field of 
the average velocity of the stellar disk under such condi- 
tions has not been calculated so far. In a general case, the 
structure of the field of the average velocity in the stellar 



systems differs substantially from the geometry of individ- 
ual trajectories of stars. Contopoulos (1978) was the first 
to calculate individual trajectories of stars in an external 
gravitational potential of spiral arms. He showed that ev- 
ery star participates simultaneously in two motions: high- 
frequency (around the center of the "epicycle") and low- 
frequency (around the Lagrangian points L2 and L4). 
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Vortex structure of NGC 157 
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Fig. 25. Restored velocity field of NGC 157 in the reference frame rotating with the pattern speed. Overlaid squares 
for each galactocentric radius show azimuthal positions of maximum value of the second Fourier harmonics of the Ha 
brightness map of the galaxy. 



As a consequence of the specific choice of the corota- 
tion radius (see the end of the previous Section), the phase 
difference between the azimuthal and radial residual ve- 
locities becomes zero at corotation, and both velocities 
simultaneously drop to zero at the centers of the vortices 
directly on the corotation circle. But, as was noted earlier, 
the qualitative picture of the restored velocity field does 
not depend on the particular choice of parameters within 



the range of their uncertainties. In all cases two anticy- 
clones appear between the spiral arms, with their centers 
located near the corotation radius. 

It should be noted that in reality the general picture 
of gas velocities in a galaxy includes not only the motions 
in the plane of a disk, but also vertical motions which con- 
tribute to the second Fourier harmonic of the line-of-sight 
velocity field (Eqs. ^l|) , and hence can be determined in 
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the frame of the model described above (Eqs.|2^, ^). The 
most convenient way to study vertical motions is to mea- 
sure gas velocities in galaxies seen nearly face-on, where 
velocity components in the plane of the disk give a small 
contribution to the line-of-sight velocity. As an example, 
we refer to the analysis of the velocity field (both from 
Hi and Ha observations) of the spiral galaxy NGC 3631, 
where vertical gas motions related to the spiral arms were 
really found (Fridman et al. 1998, 2000). 

Comparison of the location of the spiral arms with 
that of the anticyclones shows that the centers of the 
anticyclones lie between the spiral arms. Thus, the self- 
gravitational forces in the gas dominate over the forces 
of hydrodynamic pressure in this galaxy (Lyakhovich et 
al. 1996). Therefore the origin of the spiral- vortex struc- 
ture of NGC 157 must be explained in the frame of the 
gravitation concept of the density waves. 

An analysis similar to the one undertaken in the 
present work is possible only for galaxies where the nu- 
merous estimates of the observed gas velocities are well 
distributed all over the disk and spiral density waves are 
present. It demonstrates a great potential concerning the 
investigation of Fourier components of the line-of-sight ve- 
locity azimuthal distributions when a sufficient amount of 
observational data is available. 

9. Conclusions 

— The line-of-sight velocity field obtained for the spiral 
galaxy NGC 157 from interfcromctric observations in 
the Ha emission line, allowed us to obtain reliable esti- 
mates of the parameters of the disk orientation, of the 
coordinates of the dynamical center and of the rotation 
curve in the approximation of pure circular motion. 

— Deviations of the observed line-of-sight velocities from 
pure circular motion have a systematic character and 
can be described by the first three harmonics of the 
Fourier expansion of the azimuthal distribution of ve- 
locities. 

— When perturbed gas motions are taken into account, 
these do not substantially change the estimates of the 
disk inclination, the dynamical center position, and the 
major axis orientation. 

— The good correspondence between two phase curves - 
that of the modified third harmonic of the line-of-sight 
velocity field and that of the second Fourier harmonic 
of the perturbed surface density of NGC 157, gives 
evidence for the wave nature of the spiral structure in 
this galaxy. 

— The predominance of the first three Fourier harmonics 
in the line-of-sight velocity field is naturally explained 
by the prevalence of perturbations with mode m =2 in 
the density wave. This circumstance opens the way to 
restore the full vector velocity field of the ionized gas. 

— Two methods of vector velocity field restoration in a 
galaxy with a two-armed spiral are applied to the line- 
of-sight velocity field of NGC 157. They give qualita- 
tively and quantitatively similar results. Typical values 



of perturbed velocities in the density wave are 20-30 
km s~^. The maximum amplitude of the radial veloc- 
ity in the density wave is ~ 30 km s^^, that of the 
azimuthal velocity is ~ 40 km s~^. 

— Using the approximation, that the pitch angle of the 
spiral pattern is small, we propose two methods to de- 
termine the position of the corotation radius. The first 
method is based solely on the Fourier analysis of the 
observed velocity field, and gives Vc — 42" ± 5". The 
second method involves also the additional data on ra- 
dial variations of the perturbed surface density phase. 
This method gives Tc = 42" ± 6", which shows that 
these two independent approaches are in good agree- 
ment with each other. 

— The determination of the corotation radius of the spiral 
pattern and the restoration of the gas vector velocity 
field in the galaxy have enabled us to detect anticy- 
clonic structures, the existence of which had been ear- 
lier predicted analitically and from laboratory exper- 
iments on shallow water (Nezlin et al. 1986). For the 
galaxy NGC 157 we found that the anticyclone cen- 
ters are located between the spiral arms; this gives ev- 
idence for the dominance of self-gravitation over hydro- 
dynamic pressure forces near corotation (Lyakhovich 
et al. 1996). 
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Appendix A: Derivation of basic relations 

This Appendix contains a short derivation of the basic 
relations used in the Paper to close the system ( |20| ) - ( psj ) , 
in order to restore the vector velociy field of the gas, and 
to determine the corotation radius of the spiral structure. 
Whereas the final results were obtained earlier (L97; F97) 
the method of derivation presented here is new and more 
concise. 

Following the traditional way of representating small 
amplitude perturbations in the disk (see e.g. Lin & Shu 
1964; Fridman & Polyachenko 1984) we shall derive this 
relation below for the case of the tightly wound spiral ap- 
proximation. The perturbed functions for a grand design 
spiral galaxy can be presented in the following form: 

cr(r, (p, t) = ao{r) + ai{r, (f, t) = 



= ao{r) + aexpi[Jk{r)dr + 2ip-Lut-F„], (Al) 



Vrir, (f, t) = Vri{r, (/?, t) 



= Vrexpi[ J k{r)dr + 2ip-ujt ~Fr], (A2) 
V^{r, ^, t) ^ "KpoW + y^i(r, t) = 

= V^o{r) + %expi[jk{r)dr + 2ip~ujt-F^],iA3) 



where k{r) is a local radial wave number, and 



fc(r)dr > 2. 



(A4) 



The latter inequality is equivalent to the tightly wound 
arms approximation, which is in agreement with the form 
of the arms of NGC 157. 

Unlike the expressions (|^) - (|^), here we use the com- 
plex representation of the perturbed values and assume 



that the amplitudes marked by a tilde in Eqs. ( |Al| ) - (A3) 
slightly depend on r. 



Substituting (Al) - (A2) into the linearized continuity 
equation for the surface density^ 

^ + - ^(raoVri) + --^{aoV^i + aiV^o) = (A5) 
ot r or r dip 



and using (A4), we obtain 



cj(7 exp(— iFg.) = (Tofc,.T4- Gxp(— ii^r) ) 
where 

Qj = Lo — 2r2o(r') . 



(A6) 



(A7) 



In most galaxies f^o is a decreasing function of r. Thus 
according to (AT) one obtains: a) < inside corotation 
radius (at r < Vc) and a) > outside corotation radius (at 
r > Tc). Then for a trailing spiral, fcr > 0, we find from 
(Rfl) the following relation between phases: 



F„-Fr 



■K aX r < Tc 
at r > 7v 



[1 — sgn(cD)] 7r/2 . 



(AS) 



From the linearized equation of motion in the disk 
plane under the condition (A4) we obtain]] (Fridman & 
Khoruzhii, 1999): 



(A9) 



Vri Vr exp(— ii^r) . 2^0'^ 

~ % exp(-ii^^) ~ k2 ' 

® The continuity equation for the surface density, which is 
derived by integrating over z the initial 3D continuity equation, 
contains the residual velocities V\ = J povidz, which are 
residual velocities averaged over z. 

^ For barotropic perturbations the function V^o depends on 
r only (the Poincare theorem). Therefore the initial equation 
of motion can be integrated over z, to obtain the equation for 
the residual velocity V averaged over z (L97). 
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whence: 

exp[— 



^^exp[isgnMV2].(A10) 



Here k, is an epicyclic frequency, = 2ilo(2r2o + '"^^o)- As 
follows from Fig. |[ > 0, therefore from ( A10| ) we find 



Fr 



F^ = 



7r/2 at r < Tc 
-7r/2 at r > Tc 



sgn(cD)7r/2. (All) 



Adding ( |A8| ) to ( |All| ) and multiplying by -1 we arrive 
at the following relations: 



F„ = 



7r/2 at r < 
7r/2 at r > Tc 



7r/2. 



On the other hand, from Eqs. 
account (All), one can derive 



(A12) 
51) taking into 



^obs 
"3 



, sgn (cosF^) = sgn 



^obs 



(A13) 



Comparing the latter with the expressions derived 
from (|l7l), (|l|): 



Lob; 
"3 



tani^3 = sgn (cosFg) = sgn (cg 



(A14) 



we see that outside corotation the phase of the azimuthal 
perturbed velocity {F^p) coincides with the phase of the 
third Fourier harmonic {F^). 

Inside corotation the relation between phases of the 
perturbed azimuthal velocity and the third Fourier har- 
monic of the observed line-of-sight velocity depends on 
the ratio of amplitudes and C^. 

For a two-armed tightly wound spiral the ratio of am- 
plitudes of azimuthal and radial residual velocities is (see 



Eq. (|A9D) 



\Vr\ 



2ni 



{2n + rn') 
2\n-a 



ph\ 



(A15) 



Near corotation, this ratio is high because the denomina- 
tor is small. The analysis of the rotation curve of NGC 157 
(Fig. H) shows that in the inner region of the disk the value 
rd^l/dr is small, therefore 



c^/Cr ~ n/\n-n 



ph\ 



> 1. 



(A16) 



Consequently, in most of the galactic disk the phase 
of the residual azimuthal velocity will be approximately 
equal to the phase of the third Fourier harmonic of the 
observed line-of-sight velocity: 



F„ 



F. 



(A17) 



Combining the relations ( A12| ) and ( A17) we finally ob- 
tain the connection between the phases of the two observ- 
able values — perturbed surface density and third Fourier 
harmonic of the line-of-sight velocity: 



F. 



F„ 



7r/2. 



(A18) 



This relation should be valid if the grand design spiral has 
a wave nature. 

The derived relations between the phases of the differ- 
ent parameters of the perturbations can be used to close 
( pO| ) - ( p5| ) and determine the characteristics of the full 
vector velocity field of the galaxy. 



The system of equations (20) - (25) can be rewritten 



as follows: 






A3 = af'coti. 




(A19) 


B, = bf^'coti, 




(A20) 


Ar = -bf' + 6f^ 




(A21) 


B^ = + 




(A22) 


VJ^^^^+A^ = af'' + 


obs 
"3 ' 


(A23) 


VZ°'^^ + Br = af' - 


obs 
^3 > 


(A24) 


where we have inserted the sine and cosine of the resid- 
ual velocity components related to their amplitudes and 
phases: 


= d cos Fi , 




(A25) 


Bi — Ci sin Fi . 




(A26) 



The former parameter characterises the residual velocity 
variations along the dynamical major axis of the galaxy, 
and the latter describes variations along the minor axis. 



It follows from ( A19 ) - ( A22 ) that by using only ob- 
servational data, without any additional conditions, one 
can determine the vertical velocities in the galaxy to- 
gether with the components Ar and B^, which describe 
the amplitude of the perturbed radial velocity at the ma- 
jor axis of the galaxy and the amplitude of the perturbed 
azimuthal velocity at the minor axis, respectively. The ver- 
tical motions are described by the second harmonic and 
differ from a2^^ cos{2lp) + b2^^ s'm{2ip) only by a factor of 
cot? which is approximately equal to 1.3 for NGC 157. 

Two remaining equations, ( A23| ) and (A24), contain 
three unknowns, V^"'^^, A^, and Br, which cannot be 
determined without additional conditions. 



The first way to solve Eqs.(Al9) - (A24) is to use the 
relation (All) between the phases of the radial and tan- 
gential residual velocities. By using this condition, from 



( [A2l| ), ( |A22D , ( |A25D and ( |A26D we obtain 

^obs 



± 



6f 



Lobs 
°3 
Lobs 

"3 



and from ( |A23| ), (|A24| ) 



± 



^obs 



mod 2 



^obs 



-I- a' 



obs 



mod2 



(A27) 



(A28) 



The upper sign corresponds to the region inside corota- 
tion, and the lower sign - outside it. 
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By equalising the right hand terms of (A27) and (A28) 
we find at r < Tc, as well as at r > Tq'- 



''rot 



d2 



a 



obs 



and finally from (A23) and (A24) we get 

^obs (hphs luohs 
"3 \"l I "3 



Br 



(6f 76§ 



1), 



(A29) 



(A30) 



Another way to close the system of Eqs.( A19) - ( A24) 
is to use the relation ( |A8D between phases of the perturbed 
surface density and the radial residual velocity. The phase 
of the perturbed density can be determined on the basis 
of analysis of an image of a galaxy, assuming that the 
locations of the maxima of the perturbed surface density 
at a given radius coincide with the locations of the maxima 
of the sur face b right ness r elated to the spiral arms. 

From (|A2l| ) and (|A25|) we have 



Cr cos Fr 



lobs 



obs 



(A31) 



and im media tely the value Br is calculated, using its def- 
inition ( A26 ) and the phase relation (|A^): 

(A32) 



Br = Cr sin Fr = (6; 



obs 
3 



6f")tani^<, 



Substituting ( |A3^ ) into ( |A2| ), we find for both r < 
and r > Tr'. 



3d2 



whence 



2a' 



obs 



obs\ 



obs\ 



(A33) 



(A34) 



In a similar way, the relations for determination of the 
corotation radius from the observational data can be de- 
rived. 

From (|A27|) , we have 



Lobs 
"3 



^obs ^ (•^obs)2 _ (T^ohs^ 



T ^ . (A35) 



^obs _^ ^obs (^obs ^ ^obs-)2 

Taking into account that the amplitudes Cr and C^p always 
have positive values, equation ( A35| ) gives 

_ \bf'%r)\ < 0, for 



bsir)] - l&f'WI > 0, for 



r <rc, 

r > Tr. 



(A36) 



These inequalities offer a possibility to determine the lo- 
cation of corotation from the radial behaviour of \b'^^\ and 
in the observed line-of-sight velocity field. 
On the other hand, from ( A3l| ) we have 

acosF^ = T CrCosF,, = bf"" - 6f^ (A37) 



I TiObs 
I "3 



Multiplying the last equality of Eq.(A37) by cosFcr, 
we find the equation 

- 6f'')cosF^ = T acosV^, (A38) 

from which follows 

(6f^(r)-6f^(r)) cosF,(r) < 0, for r < r„ 

{bf^'ir) -bf\r)) cos F„{r) > 0, for r > r^. ^ ^ 

These two inequalities allow to find the location of coro- 
tation on the basis of the surface brightness map and the 
line-of-sight velocity field. 



